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Abstract 

Using a relativistic effective Lagrangian at the liadronic level, near-threshold 
UJ and (j) meson productions in proton proton (pp) collisions, pp ppto/cj), 
are studied within the distorted wave Born approximation. Both initial and 
final state pp interactions are included. In addition to total cross section 
data, both uj and (p angular distribution data are used to constrain further 
the model parameters. For the pp — > ppuj reaction we consider two different 
possibilities: with and without the inclusion of nucleon resonances. The nu- 
cleon resonances are included in a way to be consistent with the tt~p — > ujn 
reaction. It is shown that the inclusion of nucleon resonances can describe 
the data better overall than without their inclusion. However, the SATURNE 
data in the range of excess energies Q < 31 MeV are still underestimated by 
about a factor of two. As for the pp pp(p reaction it is found that the 
presently limited available data from DISTO can be reproduced by four sets 
of values for the vector and tensor (l)NN coupling constants. Further mea- 
surements of the energy dependence of the total cross section near threshold 
energies should help to constrain better the (f)NN coupling constant. 

PACS number(s): 25.40.Ve, 24.10.Jv, 24.30.-v, 13.75.Cs 

I. INTRODUCTION 



Heavy meson production in nucleon nucleon [NN) collisions can, in principle, provide 
important information about the short-range part of the NN interaction [1]. For exam- 
ple, the A^A^ —>■ NNuj/(j) reactions at their threshold energies probe distances between the 
two colliding nucleons of about 0.2 fm [2]. The distance corresponds to the "overlapping 
region" of the two interacting nucleons, in contrast to the distance of about 0.5 fm probed 
by much lighter pion production [3]. Therefore, investigation of such heavy meson produc- 
tion reactions should ultimately provide relevant information for testing QCD-based A^A^ 
interactions. 

Recently, there has been considerable interest in the vector mesons u and 0, in connec- 
tion with the OZI rule [4] violation, and in the strangeness content of the nucleon wave 
function [5-7]. For example, the Crystal Barrel experiments at LEAR (CERN) [8] found a 
strong violation of the OZI rule in the (f)/uj production rate in antiproton-proton {pp) annihi- 
lation. Furthermore, it was also found that the to u; production ratio in the pp —>■ pp uj / (p 



1 



reactions was enhanced by about an order of magnitude relative to the OZI prediction after 
correcting for the available phase space volume [9]. These findings may be interpreted as a 
considerable admixture of the ss configuration in the nucleon wave function. 

Another item of interest in u meson production processes is the so-called "missing res- 
onances" problem, where constituent quark models predict more states than have been 
observed experimentally [10,11]. This has been attributed to the possibility that many such 
missing resonances couple either weakly or not at all to the ttN channel, but may cou- 
ple more strongly or exclusively to the uN channel. Indeed, some theoretical studies of uj 
photoproduction were made [12,13] inspired by this possibility. 

Although, so far, no baryon resonances have been observed to decay into the uN channel, 
some theoretical efforts have been made to estimate the coupling strengths of oo to the 
experimentally observed resonances [14-16]. One of the motivations for such a study is the 
possibility to account for the observed enhancement of the low-mass dilepton pair production 
in heavy ion collisions [17]. Alternatively, the downward shift of the p and uj meson masses 
(but a smaller shift for the (f) meson mass) in the nuclear medium [18-22] is also considered as 
a possible source of the observed enhancement. Indeed if the downward shift is large enough, 
ui meson is expected to form meson- nucleus bound states [23-25]. In any case a better 
understanding of the vector meson production mechanism in free space is a prerequisite to 
study such in-medium effects, and also to study the possible couplings of the lu meson to 
(missing) resonances. However, in spite of the pronounced interest in the vector mesons p, cu 
and (f), so far there exist only a limited number of theoretical studies of the near-threshold 
pp ppuj [26-30] and pp pp(f) [26,31,32] reactions. This situation also holds for the 
near-threshold pn duo/ (j) [33,34] and pd ^^Hccj [35] reactions. 

From the experimental side, apart from the old data at high excess energies [36], only 
the total cross section data from SATURNE [37] were available until recently for pp — > ppu 
in the near-threshold region with excess energies below Q = 31 MeV, where the excess 
energy Q is defined as, Q = ^/s — 'Y^pirip with ^/s and mp being the total center-of-mass 
energy and masses of the particles in the final state, respectively. There are also total cross 
section and angular distribution data at excess energy Q = 319 MeV from the DISTO 
Collaboration [9]. Recently the COSY-TOF Collaboration has measured the total cross 
section for pp — > ppcu at two excess energies, Q = 92 and 173 MeV [38]. These fill in 
partly the energy gap between the SATURNE and DISTO data, and are critical in studying 
the energy dependence of the total cross section in the extended near-threshold region. In 
addition to the total cross sections, the COSY-TOF Collaboration has also measured the 
angular distribution of the co meson produced at Q = 173 MeV. As has been pointed out 
in Ref. [27], the angular distribution plays a major role in disentangling different reaction 
mechanisms. These new data from the COSY-TOF Collaboration, together with earlier 
data [37], offer the opportunity to investigate the pp — > ppuj reaction more in detail than has 
been done previously. Thus, in the present study we focus on the near-threshold pp — > ppuj 
reaction in free space. We study this reaction by considering two different possibilities: with 
and without the inclusion of nucleon resonances [30]. The possibility of a large, off-shell 
S'ii(1535) resonance contribution has been considered recently by the Tiibingen group [29]. 

In the present work we also consider the pp ppcj) reaction. The only data available 
for this reaction near-threshold energies are the reanalyzed total cross section and angular 
distribution from the DISTO Collaboration at Q = 83 MeV [9] . These data are not sufficient 
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to provide stringent constraints on theoretical models. In studying the pp pp(j) reaction 
we do not consider the possibility of nucleon resonances, because not enough data exist to 
either draw any meaningful conclusions about their role, or to fix new parameters associated 
with the resonances. (Also note that no baryon resonances have been observed decaying into 
the (j)N channel.) 

To study the pp — > ppuj/<p reactions, we use a relativistic effective Lagrangian at the 
hadronic level, where the reaction amplitude is calculated within the distorted wave Born 
approximation, including both the initial and final state pp interactions (denoted by ISI 
and FSI, respectively). The ISI is implemented in the on-shell approximation [2,30,39,40], 
while the FSI is generated using the Bonn NN potential model [41]. The finite width of 
the a;, which is very important near threshold energies, is also included. The up FSI is 
included only via the pole diagrams (s-channel processes). Many of the cut-off parameters 
and coupling constants necessary for the present study have already been fixed from other 
reactions in previous studies [2,27,31,39]. It turns out that the pp —>■ ppuj reaction is ap- 
parently described better with the inclusion of nucleon resonances. However, in order to 
draw a definite conclusion we need more data for exclusive observables in the energy region 
above, but close to Q = 30 MeV. This is because there is no established method to remove 
the multi-pion background associated with the a;-meson width from the raw data in order 
to extract the cross sections. The finite lu width is very important for energies Q < 30 
MeV and can possibly make the extracted data strongly dependent on the model used in 
the analysis [37] . As for the pp — > pp0 reaction, we definitely need more data to constrain 
the model parameters. 

This article is organized as follows. In Section 11 the general structure of the reaction 
amplitude in the present approach is explained. In Section 111 we discuss the pp — > ppuj 
reaction without the inclusion of nucleon resonances; in Section IV we discuss the reaction 
with the inclusion of nucleon resonances. Section V treats the pp — > pp(p reaction. The 
results are discussed and summarized in Section VI. 

II. STRUCTURE OF REACTION AMPLITUDE 

In this section we describe the structure of the reaction amplitude for the NN — > 
NNV {V — uj,(p) reaction following Ref. [39], to make this article self-contained. 

In Fig. 1 we show a decomposition of the reaction amphtude for the NN NNV {V — 
uj. (f)) reactions. The reaction amplitude is calculated in the distorted wave Born approxima- 
tion using a relativistic meson exchange model. We begin by considering the mcson-nucleon 
(MN) and A^A^ interactions as the building blocks. Then, we consider all possible combina- 
tions of these building blocks in a topologically distinct way, with two nucleons in the initial 
state, and two nucleons plus a meson in the final state. Here diagrams leading to double 
counting, e.g., those contributing to mass and vertex renormalizations must be excluded, 
since we use the physical masses and coupling constants. The ellipsis in Fig. 1 indicates those 
diagrams that are more involved, or higher orders, which are not included in this work. In 
particular, we neglect the MN FSI, which otherwise would be generated by solving the 
three-body Faddeev equation. 

In order to make use of the available potential models of NN scattering, we perform 
the calculation within a three-dimensional formulation which is obtained from the Bethe- 
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Salpeter equation by maintaining the relativistic unitarity and Lorentz covariance of the 
resuhing amphtude. We foUow the procedure of Blankenbecler and Sugar [42] as adopted 
in the Bonn NN potential model [41]. The vector meson production amplitude, M, may be 
written as 

M = (1 + Tj-^hG^f^*){e* • J)(l + iCj+^r/^^) , (1) 

where T(^ij) stands for the NN T-matrix in the initial(i)/final(/) state, is the three- 

dimensional Blankenbecler-Sugar (BBS) propagator, and e* is the polarization vector of 
the vector meson produced. The superscript ± in T(^ij) and G'(jj) indicates the boundary 
conditions, (— ) for incoming and (+) for outgoing waves. The production current J^, which 
is the MN T-matrix with one of the meson legs attached to a nucleon (first diagram on the 
r.h.s. in Fig. 1) is defined by 

J'^^Yl [T{MN^M'N)] 1 iPm' [^1,'Nn]2 + (1^2), (2) 
M' 

where T(^mn^m'n) stands for the MN T-matrix describing the transition M'N — > MN, and 
^M'NN ^M' are respectively the M'NN vertex and the corresponding meson propagator. 
The subscripts 1 and 2 denote the two interacting nucleons 1 and 2. The summation is 
over the intermediate mesons M'. We note that if the four dimensional full two-nucleon 
propagator is used, the reaction amplitude given by Eq. (1) would have an additional term 
to avoid the double counting arising from the term involving iG-^^T;^ and/or iG^j^^T^ and 
the current J^, since also contains meson-exchange NN interactions. This additional 
term vanishes when wc use the reduced three dimensional propagator G{.iif) which restricts 
the energy of the propagating two nucleons to be on their respective mass shells. 

In the near-threshold energy region, the two nucleon energy in the final state / is very 
low and hence the NN FSI amplitude, T^"^^ in Eq. (1), can be calculated from a number 
of realistic A^A^ potential models. In the present work we use the Bonn NN potential 
model [41]. This model is defined by a reduced three dimensional BBS version of the Bethe- 
Salpeter equation, 

T = V + ViGT , (3) 

where G denotes the BBS two-nucleon propagator, consistent with those appearing in 
Eq. (1). (Note that the factor —i difference in the definitions of V and T from those in 
Ref. [41].) 

For the A^A^ NNuj/(f) reactions the A^A^ initial state interaction (ISI) amplitude, Tj*-"*"^ 
in Eq. (1), must be calculated at incident kinetic beam energies above 1.89 and 2.59 GeV, 
respectively. There exists no accurate A^A^ interaction model with which one can perform 
calculations reliably at such high incident beam energies. In the present work we follow 
Ref. [40], and make the on-shell approximation to evaluate the ISI contribution, which was 
also applied in the study of the A^A^ NNr) reaction [39] . This amounts to keeping only the 
5— function part of the Green function C; in evaluating the loop integral involving iG-^-'T^*-^^ 
in Eq. (1). The required on-shell A^A^ ISl amplitude is calculated from Ref. [43]. As discussed 
in Ref. [40], this is a reasonable approximation to the full A^A^ ISI. In this approximation 
the basic effect of the A^'A'" ISI is to reduce the magnitude of the meson production cross 
section [39,44]. In fact, it is easy to sec that the angle-integrated production cross section 
in each partial wave state j is reduced by a factor Xj [40] : 
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= ,7,(p) cos\5,{p)) + i[l - r^M' < i[l + r^,{p)f , (4) 

where Sj{p) and ?7j(j9) denote respectively the phase shift and corresponding inelasticity, 
and p is the relative momentum of the two nucleons in the initial state. The A^A^ ISI has 
been considered fully by Batinic et al. [44] in the study of the pp — > ppr] reaction, whose 
threshold corresponds to an incident energy of about 1.25 GeV. Their results support the 
on-shell approximation used in the present work. Quite recently, Baru et al. [45] have also 
investigated the effects of the A^A^ ISI in the A^A^ — > NNt] reaction. Although there are 
obviously (off-shell) effects which are absent in the on-shell approximation, the results of 
Ref. [45] also show that the bulk of the A^A^ ISI effect is accounted for in the on-shell 
approximation. In particular, we do not expect the off-shell effects of the ISI to change the 
conclusion of the present work. 

Next, we consider the production current J'^ defined by Eq. (2) based on meson exchange 
models. Following Refs. [31,39,46], we split the MN T-matrix in Fig. 1 into the pole {T^^) 
and non-pole (Tj^^) parts and calculate the non-pole part in the Born approximation. Then, 
the MN T-matrix can be written as [47] 

Tmn — T^j^ + > (5) 

where 

Tmn = XI fliNBWBfMNB , (6) 

B 

with J'mnb and qb denoting the dressed meson- nucleon-baryon (MNB) vertex and baryon 
propagator, respectively. The summation is over the relevant baryons B. The non-pole part 
of the T-matrix is given by 

rpNP _ yNP , yNP-nr^NP /yN 
^ MN — ^ MN ^ ^ MN''^-'- MN 1 {' ) 

where = Vmn — Vmni with V^^v denoting the pole part of the full MN potential Vmn- 
V^j^ is given by an equation analogous to Eq. (6) with the dressed vertices and propagators 
replaced by the corresponding bare vertices and propagators. We neglect the second term of 
Eq. (7) and hence the full MN T-matrix in Eq. (2) is approximated as Tmn — ^mat + ^mn- 
With the approximation described above, the resulting current consists of baryonic 
and mesonic {J mec) currents. The baryonic current is further divided into the nucleonic 
{Jnuc) nucleon resonance ( J/4s) currents, so that the total current is written as 

The vector meson production currents are illustrated diagrammatically in Fig. 2, where 

V stands for the or meson. Note that they are all Feynman diagrams and, as such, 
they include both the positive- and negative-energy propagation of the intermediate state 
particles. The nucleonic current is constructed consistently with the A^A^ potential in Eq. (3). 
For the mesonic current, it turns out that we may consider only the V pi: {V — uJi(t>) 
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exchange- current contribution for both the pp — > ppu/<f) reactions. The nuclcon resonance 
current J^^^,, inchidcd only in the pp ppuo reaction, is explained in Section IV. 

Here, some general remarks on the meson production currents are in order [39]. As a 
consequence of using a three-dimensional reduction of Bethe-Salpeter equation, the definition 
of the energy (time component) for the intermediate state particles in the production currents 
becomes ambiguous. In order to be consistent with the NN interaction used in the present 
work, we follow the BBS three-dimensional reduction prescription: (1) The energy of a 
virtual meson at the MNN vertex is taken to be go = ^(p) ~ where e{p)[e{p')] is the 

energy of the nucleon before [after] the emission of the virtual meson, with e{p) = ^/p^ -|- m^. 
(2) The energy of the intermediate state baryon B in the nucleonic and resonance currents 
is taken to be pq — uj{k) -\- s{p') at the B ^ M + N vertex, while at the N ^ M + B vertex 
it is taken to be pq = e{p') — uj{k), where uj{k) is the energy of the meson produced in the 
final state. The BBS reduction prevents three particle cuts which occur in a more exact 
calculation. 



III. PP PPu) WITHOUT RESONANCE 

In this Section, we consider the pp —>■ ppuj reaction without the inclusion of nucleon 
resonances. Then, the total cj-meson production current J'^ may be given by the sum of the 
nucleonic and ujp-K meson-exchange currents, = Jn^c'^ Jmec-i shown in Fig. 2 {V = uj). 

The nucleonic current J^^^ is defined by 

JLc = E K^^.-^ + UiS.V^) , (9) 

with denoting the cvNN vertex and Sj the nucleon (Feynman) propagator for the nucleon 
j. The summation is over the two interacting nucleons, 1 and 2. U stands for the meson- 
exchange A^A^ potential which is, in principle, identical to the driving potential V used in 
the construction of the A^A^ interaction (Eq. (3)) , except that here meson retardation effects 
are retained following the Feynman prescription. 

The CO production vertex loNN, in Eq. (9), is obtained from the Lagrangian density, 

where m^v, i^Nix) and cl-'^(x) stand for the nucleon mass, nuclcon and c<j-meson fields, re- 
spectively. QujNN denotes the vector coupling constant and K^i = fujNN /Ql^nn {dujNN 0), 
with fuiNN the tensor coupling constant. J^„^ defined by Eq. (9) is illustrated in Fig. 2a. 

As in most meson-exchange models of hadronic interactions, each hadronic vertex is 
accompanied by a form factor in order to account for the composite or finite-size nature 
of the hadrons involved. Thus, the uNN vertex obtained from the above Lagrangian is 
multiplied by a form factor [31,39], 

A^ 

^-^^(P"^ - A^^ + ^/^rnir ' ^^^^ 
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where is the four-momentum squared of either the incoming or outgoing off-sheU nucleon. 
It is normahzed to unity when the nucleon is on its mass shell, = mj^. Following Ref. [31], 
we adopt g^jNN = 9.0 in Eq. (10). The vector to tensor coupling constant ratio, k^,, is not 
well established; in fact the values quoted in literature are relatively small and vary in a 
range, — 0.16±0.01 < Ko; < +0.14±0.20 [27]. Therefore, we consider and Ajv, respectively 
in Eqs. (10) and (11), as free parameters in the present work. 

The ujpiT vertex for u production in the meson-exchange current, J^^^c (^ig- 2b), is derived 
from the Lagrangian density, 

jC^pnix) = -^^e^^,^d-pP(x) ■ d''7r{x)u;^{x) , (12) 

where Sapup. is the totally antisymmetric Levi-Civita tensor with £0123 = —1- The cufm vertex 
obtained from the above Lagrangian is multiplied by a form factor, 

where A^^ = Ap = A^ is assumed [31]. It is normalized to unity at = and = m^, 
consistent with the kinematics at which the coupling constant g^jp.,^ is extracted. 
The meson-exchange current is given by 

JLc = [^^iViv(?p)]l^^aM?p)^S^J?P, Qn, k^)tA{q^)[r^NN{q.)]2 + (1^2), (14) 

where Da/siqp) and A(57r) stand for the p- and 7r-meson (Feynman) propagators, respectively. 
The vertices F involved arc self-explanatory. The coupling constant, g^^^p^ = 10.0 in Eq. (12), 
has been fixed from a systematic study of pseudoscalar and vector meson radiative decays 
combined with the vector meson dominance assumption [31,48]. Its sign is fixed from a study 
of pion photoproduction in the 1 GeV energy region [49]. The pNN and ttNN vertices in 
Eq. (14) are consistent with those in the Bonn-B NN potential [41], except that here we 
use the pseudovcctor-coupling, which is consistent with the chiral constraints in the lowest 
order [50], instead of the pseudoscalar-coupling for the nNN vertex. In addition, the cut- 
off parameter, A^^nn = 1300 MeV, is adopted at the ttNN vertex. (One could use the 
pseudoscalar-coupling, or even the admixture of the pseudoscalar- and the pseudovector- 
couphngs [51]. Note that these two couphngs entail different orders in chiral counting [50]. 
A test calculation by the use of the pseudoscalar-coupling with the same model parameters, 
turned out to give an enhancement of the total cross section by about a factor of ten near 
the threshold, and the enhancement became far larger as the excess energy Q increases.) 
We are, then, left with the cut-off parameter Ap — A^ in Eq. (13) which will be treated as 
a free parameter in the present work. 

Next, we explain how the model parameters k^^, A^r, and Ap = A^ arc determined in the 
present approach. In Ref. [27] it was pointed out that the angular distribution of the emitted 
u mesons is a sensitive quantity for determining the absolute amount of nucleonic as well as 
mesonic current contributions in addition to the relative sign of the two amplitudes. This was 
also demonstrated in Refs. [30,52]. (The value quoted in Ref. [52] should read Qunn — +9.0.) 
Furthermore, the shape of the a; angular distribution is particularly sensitive to the value of 
Ki^, the tensor to vector coupling ratio [30,52]. Thus, we can make use of both the cu angular 
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distribution and total cross section data from COSY-TOF [38] at Q = 173 MeV to fix these 
three parameters. We obtain a reasonable fit to the data with the values, Ap = A^^ = 1000 
MeV, Ajv = 1190 MeV and K^J ~ —2.0. Table I summarizes all the parameter values within 
the approach in this Section. Here, it should be mentioned that, aXQ ~ 173 MeV, the energy 
involved in the final pp subsystem extends beyond the pion threshold. The NN FSI used in 
the present work has been developed to fit the phase-shifts only up to the pion production 
threshold. Therefore, strictly, we are beyond the applicability of this interaction. However, 
the final pp energy involved is not large enough to introduce any significant deviation. 

The dependence of the u angular distribution is illustrated in Fig. 3. Here, since the 
mesonic current contribution gives a flat angular distribution, and because we are interested 
in the k^. dependence of the shape, we have kept the mesonic current contribution unchanged 
(Ap = Att = 1000 MeV) and varied both and K^. The latter parameter has to be 
varied in order to keep the total cross section fixed. The results show indeed the shape is 
sensitive to the value of k^. Note in particular that values ol k^^ > —1 are clearly unable to 
reproduce the data. We note here that the value of the tensor coupling, f^^NN — 0, in the 
Bonn NN potential used in the FSI in the present work, is quite different from the value 
of fuiNN = i^LjdujNN — —2.0 X 9 = —18 used to reproduce the angular distribution data. 
However, the exchanged lv meson in the Bonn A^A^ potential model is associated with the 
virtual lu meson and represents the isoscalar- vector quantum numbers exchanged, and is not 
necessarily related to the physical cu meson. 

Next, in order to see how the shape of the u angular distribution is sensitive to the value 
of gcjNN, we keep the value, k^j = —2.0 fixed, which reproduces the u angular distribution 
very well, and calculate the u angular distribution for different values of QujNN- Here again 
the mesonic current contribution has been kept fixed. We show in Fig. 4 (the right panel) 
the result obtained using the value, {QojnnY /^t^ — (17.37)^/47r = 24, approximately the 
value used in the Bonn NN potential, together with one of the reasonable fits (the left 
panel) obtained with gu)NN = 9.0. Recall that the total cross section is normalized in 
both calculations. The result obtained with {q^nnY / '^'^ = 24 also gives a good fit to the 
data, where the change in the coupling constant Qojnn, 9.0 — > 17.37, is compensated by 
the change in the cut-off parameter Ajv, 1190 1020 MeV. This implies that, for a given 
mesonic current contribution, the shape of the u angular distribution is not sensitive to the 
value of QcoNN, but is sensitive to the value of k^^. This implies that the different momentum 
dependence introduced via the tensor coupling to the nucleonic current plays an important 
role to the shape of the uj angular distribution. 

With all the parameters fixed, we next study the energy dependence of the total cross 
section. In Fig. 5 we show the predicted energy dependence of the total cross section 
with various effects: (1) effects of the finite u width denoted by "width", (2) the initial 
and final state pp interactions denoted by "ISI" and "FSI", respectively, and (3) using a 
constant matrix element denoted by "Phase space", where "(e* ■ J) = constant" is used in 
Eq. (1). The results with "ISI-^FSI-Fwidth" in Fig. 5 (the bottom-right panel) include all 
the effects considered in the present section, and should be compared with the data. Recall 
that the present parameters are adjusted so as to reproduce both the total cross section 
of SO.Sfib and the u angular distribution at Q = 173 MeV. In Fig. 5, we also show the 
result of Ref. [37] denoted by "Hibou et al." used in the analysis of the SATURNE data. 
Obviously, the present result underestimates the SATURNE data [37] to a large extent in 
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the range of excess energies, Q < 31 MeV. The reason for this discrepancy can be attributed 
to an overestimate of the mesonic current contribution as fixed at Q = 173 MeV. This 
resuhs in a substantial reduction of the cross section close to threshold, due to a much 
stronger destructive interference between the nucleonic and mesonic current contributions 
as the excess energy decreases. In fact, the large overestimation of the mesonic current can 
be verified from the 7r~p — > oon reaction. With all the parameters fixed to reproduce the 
pp ppuj at Q = 173 MeV, we have looked at the model prediction of the energy dependence 
of the 7r~p un total cross section. It turns out that the model largely overestimates the 
data from Ref. [53] as the center-of-mass energy W increases due to the rapidly increasing 
cupTT exchange current contribution. We were unable to reproduce both the pp — > ppcu data 
at Q = 173 MeV and the energy dependence of the 7r~p cun total cross section within our 
approach which considers only the nucleonic and mesonic currents. 

In an attempt to improve the agreement in the present approach, we have investigated 
the effect of different form factors for the vector and tensor couplings at the uj production 
vertex in the nucleonic current. Different momentum dependences of the vector (Fi) and the 
tensor (F2) strong form factors at the uNN meson production vertex are quite possible. In 
particular, recent experimental results from the Jefferson Lab [54] for the ratio of the proton 
electric and magnetic form factors, iJ^pG Ep{Q^) / G MpiQ"^) {l^p- proton magnetic moment), 
show a linear decrease as the four-momentum transfer squared (Q^) increases, namely, the 
vector (Fi) to the tensor (F2) electromagnetic ratio is, F2{Q'^) / Fi{Q^) ~ 1/Q [55], which 
shows that Fi and F2 have different momentum dependences. (The perturbative QCD 
predicts F2{Q'^) / Fi^Q"^) ~ l/Q^ [56].) Using different cut-off values for the vector and 
tensor form factors, respectively, Aat^ = 1300 MeV and Aat^ = 1600 MeV, and using the 
same functional form of Eq. (11), we have recalculated the energy dependence of the total 
cross section. Note that these parameters together with other parameters, Ap = = 1450 
MeV and n^i = —0.5, can reproduce the uj angular distribution data ai Q = 173 MeV 
reasonably well. However, the predicted cross section is enhanced by only 10 ~ 20 % at near 
threshold energies, and still underestimates substantially all the SATURNE data points [37]. 
We have also considered a possible contribution of the ujaa mesonic current by assigning a 
reasonable range of values for the coupling constants and cut-off parameters associated with 
this current. But this also gives only a small contribution. Thus, within the approach of 
considering only the contributions from the nucleonic and mesonic currents, it seems unlikely 
to be able to reproduce the measured energy dependence of the total cross section in the 
range of excess energies Q < 173 MeV. Of course, we could fix the model parameters by 
fitting the total cross section at a lower excess energy point. However, the model would then 
overestimate the total cross sections near Q — 173 MeV substantially, and also it would be 
difficult to reproduce the uj angular distribution at Q = 173 MeV. 

Apart from the difficulties mentioned above, we also note that the rather large (negative) 
value of K.i_j = —2.0 required to reproduce the lj angular distribution, is not easily reconciled 
with other nuclear processes. For example, such a large value of leads to a rather strong 
NN (isoscalar) tensor force. This will affect the NN tensor force, given primarily by the 
TT and p meson exchange, in such a way that it becomes extremely difficult to describe 
NN scattering and deuteron properties. In fact, a rough calculation [57] shows that it is 
nearly impossible to describe the NN phase-shift data with such a strong tensor coupling 

{fu>NN — K'u>9ojNN — —18). 
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Thus, although there is currently no definite experimental evidence for the lo meson to 
couple to any nucleon resonance, it would be natural to expect that some resonance currents 
give contributions, since high nucleon incident energies are involved in the near-threshold 
Lo production in pp collisions. The reduction of the mesonic current at high excess energies 
may then be compensated by the nucleon resonance current contributions. We study such 
a possibility in Section IV. 

IV. PP PPuj WITH RESONANCE 

In order to limit the number of resonances considered and thereby avoid the introduction 
of an excessive number of new parameters, we restricted the resonances to those: (1) that 
appreciably decay to the A'"-|-7 channel so that the vector meson dominance (VMD) assump- 
tion may be used to produce a;, (2) whose mass distributions confined around (mAr-t-m^,) and 
therefore contribute maximally at near threshold energies, (3) that can describe consistently 
the iT~p uon reaction. In addition, to see if a dominant 5*11(1535) resonance contribution 
as reported in Ref. [29] is possible, and since many parameters associated with 5*11(1535) are 
under better control than those for other higher resonances, we also include this resonance 
in the present study. As a result, we consider contributions from the following four nucleon 
resonances, 5(1535)(i") * * * *, P(1710)(i^) * **, D(1700)(|^) * ** and P(1720)(|^) * * * *, 
where we list the spin, parity and status of the corresponding resonances explicitly [58]. The 
construction of the resonance current, and the associated details are given in Ref. [39]. 

The resonance current J^^g contribution to the pp — > ppuj reaction arises from the spin- 1/2 
('^V2res) spin-3/2 {Jy2res) rcsonancc currents in the present approach: 

Jres = Jl/2res + '^3/2res ' (^^) 

The spin- 1/2 resonance current, in analogy to the nucleonic current, is defined by 

Jl/2res = E E {KjNAMr,^ + t/iv*^-S^^*^^^.) . (16) 
i=l,2 N* 

Here r[_^^-^* stands for the ujNN* vertex involving the nucleon j. Sn*{p) = {'^ + fnN*) / {p^ — 
m%t +imN*TN*) is the A^* resonance propagator, with tun* and Tn* denoting the mass and 
width of the resonance, respectively. The summation is over the two interacting nucleons, 
j = 1 and 2, and also over the spin-1/2 resonances, N* = 5*11(1535) and Pii(1710). In 
Eq. (16) Un* {Un*) stands for the A^A^ NN* {NN* NN) meson-exchange transition 
potential, and is given by 

Un* ^ E ^MJVA^*(g)^AM(g^)^MiViv(g) + E ^MNN*(Q)^Di^^^M)(q)^MNN((l), (17) 

M=7r,i7 for Sii (1536) M=p,U) 
M=7r,?),CT for P]^]^ (1710) 

where Am(q'^) and Dh^(^m){(i) arc the (Feynman) propagators of the exchanged pseudoscalar 
(scalar) and vector mesons, respectively. Tmnn{(i) and r^^^(g) denote the pseudoscalar 
(scalar) and vector MNN vertex, respectively. These vertices are taken consistently with 
the A^A^ potential V appearing in Eq. (3), except for the type of coupling at the nNN vertex 
and the cuNN coupling constant. (Recall that we use the pseudovector-coupling instead of 
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the pseudoscalar-coupling at the nNN vertex.) The exception also apphes to the spin- 
3/2 resonance current, Jy2res- analogous expression to Eq. (17) also holds for Un*. In 
Eq. (17) we have an extra, cr-meson exchange in the A^Pii(1710) ^ NN transition potential, 
since the Pii(1710) Ntitt decay branch is relatively large [58]. Thus, we simulate these 
two pions conveniently by a cr-meson. 

Following Refs. [14,39,59,60], the transition vertices Tmnn* and I'mnn* i^^) 
and (17) for spin-1/2 resonances are obtained from the interaction Lagrangian densities. 



zAr(±) + 



1- A 



niN* ± TTlN J ^ 

1- A ^ 



r}{x) \ iPn{x) + h.c. , (18a) 



J^aNPiA^) = gaNPii'4>Pii{x)aijjN{x) + h.C. , 



T ■ 7r{x) } ijjNix) + h.c. , (18b) 

(18c) 



■'pNN- 



^x) 



rriN* + niN 
+ h.c. , 



rriN* + TTiN 
+ h.c. , 



i)N*{x) 



IpN* {x) 



p(=F)a2 

^i- + rW i-td, + n^a.^d'') 

niN* + rriN 



rriN* + TTi^ 



a;^(a;) ) i)N{x) 

(18d) 

f-p^'ix) \ijN{x) 

(18e) 



where 7r(a;), a;^(a;), p^{x) and ipN*{x) denote the tt, o;, p and spin-1/2 nucleon resonance 
fields, respectively. The upper and lower signs refer to the even(-l-) and odd(-) parity reso- 
nances, respectively. The operators F^^^ and F^/f'' in Eqs. (18a) - (18e) are defined by 



(19) 



The parameter A in Eqs. (18a) and (18b) controls the admixture of the two types of couphngs: 
pseudoscalar (ps-coupling: A = 1) and pseudovector (pv-coupling: A = 0) for an even parity 
resonance and, scalar (A = 1) and vector (A = 0) for an odd parity resonance, where both 
choices of the parameter A give equivalent results when baryons are on their mass shells. 
In this work we take A = 0. Note that in principle we should not allow only the pure 
r^''^^ couphng in Eqs. (18d) and (18e), because unlike the VNN vertex (V ^vector meson), 
this coupling alone at the VNN* vertex prevents us from estimating its strength using the 
VMD, since it violates gauge invariance. In the present work, we use a more general gauge 
invariant Lagrangian density as used in Ref. [60] based on Ref. [14]. 

Similar to the case of spin-1/2 resonances, the spin-3/2 resonance current is defined by 



7M 

'^3/2res 



(20) 



3=1,2 N* 



Here T^^jy, stands for the cjNN* vertex function involving the nucleon j. Saf3{N-'){p) = 

+ niN*) {-gap + 7a7/3/3 + {laPp " Palfi) /^rUN* + 2paPf3/3m%*} / {p'^ - 171%* + imN*FN*) 

is the spin-3/2 Rarita-Schwinger propagator. The summation is over the two interacting 
nucleons, j — 1 and 2, and also over the spin-3/2 resonances, N* — Di3(1700) and Pi3(1720). 
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In Eq. (20) U^, (U^,) stands for the NN NN* {NN* NN) meson-exchange transition 
potential, and is given by 

{q) + E ^MNN* 

^mnn{q) ) (21) 

M=7r,?j M=p,u> 

where ^MNN*iq) ^'MNN'i^) denote the pseudoscalar and vector MNN* vertices, respec- 
tively. An analogous expression to Eq. (21) also holds for U^». 

The MNN* vertices involving spin-3/2 nucleon resonances in Eqs. (20) and (21) are 
obtained from the Lagrangian densities [39,59], 

^%*{^) = (^-^) rN*{^)Q,,{z)T(^^i^r,{x)d^v{x) + h.c. , (22a) 
4^iv*(^) = (^-^) rN4x)Q,.{z)T^^^fM^) ■ d'nix) + h.c. , (22b) 
>cSiv*(^) = T^ (^H^j ^j^.(a;)e,.(;.)r?Viv(x)a;^^(a;) 

- (^7^ j {dxrN^{x)Q,.{z)T(^^M^)) uj^'^ix) + h.c. , (22c) 
4Siv.(^) = j rN*(x)e,^{z)rt^rM^) . p>^'^(x) 

' (2) \ 

^ ] {dy^V^.{x)Q,MV^^^Ti,^{x)) ■ p^^{x) + h.c. , (22d) 

where e^^(z) = g^,^ - (^ + l/2)7;,7,., and u^''{x) = d^uj^ix) - d'^u^x) and p^^(x) = 
d^p^{x) — d'^p^{x). In order to reduce the number of parameters, we take z = — 1/2 in the 
present work. 

Following Ref. [39], the relevant coupling constants associated with the resonance cur- 
rents are calculated utilizing the Particle data [58] whenever available; they are determined 
from the centroid values of the extracted partial decay widths (and masses) of the resonances. 
Those couplings involving vector mesons, are estimated from the corresponding measured 
radiative decay width in conjunction with the VMD, although uncertainties in the data are 
large. In order to reduce the number of free parameters, the ratio of the VNN* {V = p, u) 
coupling constants for the spin-3/2 resonances, £'13(1 700) and Pi3(1720), have been fixed 

to be 9vNN*/9vNN* = -2-1, the same ratio as that for ^iVPgs (1232) Ma^Pss (1232) = "2.1, 
extracted from the ratio of E2/M1 = —2.5% determined from pion photoproduction mea- 
surements [61]. 

Following Refs. [39,46], and in complete analogy to the nucleonic current, we introduce 
the off-shell form factors at each vertex involved in resonance currents. We adopt the same 
form factor as given by Eq. (11), with mN replaced by tun* at the MNN* vertex, in order 
to account for the N* resonance being off-shell. The MNN* vertex, where the exchanged 
meson is also off-shell, is multiplied by an extra form factor Fm{q^) in order to account for 
the meson being off-shell (see Eqs. (17) and (21)). The corresponding full form factor is. 
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therefore, given by the product FN{p'^)FM{q'^), where M stands for the exchanged meson 
between the two interacting nucleons. The form factor FM^q^) is taken consistently with the 
NN potential V in Eq. (3); the only two differences are the normalization point of Fp^^^{q^) 
and the cutoff parameter value of F^{q^). Here, the form factor for vector mesons Fp^^j{q^) 
is normalized to unity at = in accordance with the kinematics at which the coupling 
constant Qp^ujNN* was extracted, i.e., = (A^^^/(Ap^j^ — q'^))^. For the pion form factor 

FT^{q'^) we use the cutoff vahie A^ = 900 MeV. To be consistent in this section, this value 
(rather than A^ = 1300 MeV) is also used in the form factor at the ttNN vertex appearing 
in the mesonic current constructed. 

Since we now include the resonance current, the free parameters in the nucleonic and 
mesonic currents in the previous Section have to be readjusted. In addition, we will consider 
the (jNN* coupling for N* = Pii(1710) as a free parameter in the present work. In order 
to fix these free parameters in this Section, we use the 7r~p — > um total cross section data 
from Ref. [53] (see Ref. [15] for a discussion about the data), in addition to the pp — > ppcu 
total cross section and angular distribution data from the COSY-TOF Collaboration [38]. 
We note that at the excess energy of Q = 173 MeV for the pp — > ppu reaction, the center-of- 
mass energy W of the subsystem 7r~p — >■ un appearing as a building block in the description 
of the pp — > ppuu reaction, will reach a maximum value of ~ 1.9 GeV. Thus, we fix the 
parameters so as to reproduce the measured energy dependence of the Tr~p um total cross 
section data up to ~ 1.9 GeV. 

We show in Fig. 6 the calculated energy dependence of the total cross section obtained 
with a selected parameter set. At lower energies W , Z^i3(1700) and Pi3(1720) contributions 
are dominant, but neither S'ii(1535) nor Pii(1710) give appreciable contributions for the 
7r~p — >• uin total cross section. Furthermore, many trial calculations show that, without 
including the resonances it is very difficult to reproduce the near-threshold behavior of 
the 7r~p — > (jjn total cross section up to energies 1^ ~ 1.9 GeV using a reasonable set 
of parameters. However, with the inclusion of the resonances, we need a softer (stronger) 
form factor for the upn vertex to fit the overall energy dependence exhibited by the data. 
In particular, the part of the form factor which accounts for the off-shell behavior of the 
exchanged p meson requires a dipole form, 

F^^{ql ql) ^ FMl)FM) = (^5^) (f^) ' ^^3) 

with Hp = 2,Ap = 850 MeV and A„ = 1450 MeV. A cutoff parameter vahie of = 1100 
MeV has been also determined at the uNN meson production vertex. A different form 
factor, exp{(3qp) exp(— al^^), was introduced at the ujpn vertex in Ref. [62] to overcome the 
difficulties in reproducing the data using the monopole form factor Fp{q'^) in Eq. (23). 

Next, using the lo angular distribution data from COSY-TOF [38], we further fix pa- 
rameters associated with the Pii(1710) resonance, namely, the coupling constant QaNPn 
associated with the a exchange introduced effectively to simulate the observed decay chan- 
nel, Pii(1710) N+2tt. The value for the coupling constant QcrNPn is adjusted to reproduce 
the pp — > ppu total cross section of 30.8/i6 a.t Q = 173 MeV. Thus, in this procedure, the 
value obtained for QaNPn is not strictly related to the branching ratio for the N-\-2t: channel; 
instead, the contribution from Pii(1710) should be regarded as also taking into account the 
other possible resonance contributions not included explicitly in our model. 
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We show in Fig. 7 the u angular distribution calculated by fitting the coupling constant 
9aNPii to the total cross section of 30.8/^6, together with a more reasonable value of = 
—0.5. Recall that with the value of ~ —2.0 obtained in Section III, it would be very 
difficult to describe the NN scattering data consistently [57]. Two values for QaNPu are 
found to be able to reproduce the total cross section of 30.8/^6 at Q = 173 MeV. The results 
for the u! angular distribution are shown in Fig. 7, for those obtained with gaNPn — —4.3 
(the upper panel) and QaNPn = +4.8 (the lower panel). Although the value. QaNPn = +4.8 
reproduces the uj angular distribution data from COSY-TOF [38] better, the result for the 
energy dependence of the total cross section is worse than that with QaNPn — —4.3. Thus, 
we will show only the results obtained with QaNPu — —4.3. We summarize in Table II all the 
parameters fixed in the present approach, i.e., with the inclusion of the nucleon resonances. 

Next, in Fig. 8 we show the energy dependence of the pp ppuo total cross section 
calculated using the fixed parameters in Table II. The result is greatly improved compared 
to that without the inclusion of any nucleon resonances studied in Section III. (See Fig. 5.) 
However, it still underestimates the SATURNE data [37], which are in the range of excess 
energies Q < 31 MeV, by about a factor of two. Thus, further investigation is needed to 
understand better the near-threshold pp ppuj reaction. As already mentioned, we also 
need more data for exclusive observables in the energy region above but close to Q = 30 MeV, 
because there is no established method for removing the multi-pion background associated 
with the a;-meson width from the raw data to extract the cross sections. The effect of the 
width is very important in the energy region, Q < 30 MeV, and the extraction of the cross 
section can become highly model dependent. 

In order to explore the sensitivity of more exclusive observables than cross sections to the 
nucleon resonances in the pp — > ppuj reaction, we have also calculated the spin correlation 
functions and analyzing power at Q = 92 and 173 MeV, with and without the inclusion 
of the nucleon resonance currents. Here, we just mention that although some of the spin 
correlation functions exhibit some sensitivity to the presence of nucleon resonances, judging 
from the CTirrcntly achieved precisions for the pp ppuj data, such a presence may not 
be sensitive enough to be distinguished experimentally. A more thorough and complete 
study of the role of nucleon resonances in the production of ui mesons in NN collisions, 
especially in a combined analysis of the pp — > ppu and pn du reactions, will be reported 
elsewhere [63]. Such an analysis will consider, not only the spin observables, but also the 
invariant mass distributions. The pn duj process, for which the total cross section data 
have been reported quite recently [64], will provide additional constraints on the model 
parameters. 

Considering the results shown in Fig. 8, one possibility for improving the agreement with 
the data would be to introduce extra resonances, which enhance the total cross section at near 
threshold energies but only moderately enhance at excess energies around Q = 173 MeV, 
if such adequate candidates exist. On the other hand, the introduction of new resonances 
would introduce more ambiguities. The other effect to be investigated is the loN FSI, which 
is expected to enhance the total cross sections at near threshold energies, because a QCD 
sum rule study of the meson-nucleon spin-isospin averaged scattering lengths for the vector 
mesons p,uj and 0, suggests attractive VN (V = p,u,(f)) interactions [21]. 

Next, in Fig. 9 we show a decomposition of each resonance contribution to the energy 
dependence of the pp — > ppcu total cross section. Close to threshold energies the dominant 
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contribution comes from Di3(1700), while at higher excess energies, the dominant contri- 
bution comes from Pii(1710). although the contribution of this resonance was negligible in 
the 7r~p — > ujn reaction. Here, again the 5*11(1535) resonance contribution is very small in 
our model. We should mention that we also studied the contribution from the 5'ii(1650) 
resonance in the present approach, but it did not give an appreciable contribution. Thus, 
if we want to be consistent with both the 7r~p — > oun and pp ppuo reactions, it appears 
necessary to include at least three nucleon resonances, Pii(1710), ^13(1700) and Pi3(1720), 
in the present approach. 

Finally, before leaving this section, we should mention that the treatment of the 7r~p 
ujn reaction should be improved in our model. In particular, as shown by Penner and 
Mosel [15]. effects of higher-order terms than the Born term in the t-matrix equation is 
important and, consequently, they should be taken into account in a better way than through 
a form factor as has been done in the present work. 

V. PP -> PP(p WITHOUT RESONANCE 

The pp — > ppcf) reaction can be treated in an analogous way to the pp — > ppcu reaction 
in Section III, namely, considering only the nucleonic (Jnuc) mesonic {J^ec) current 
contributions. However, the scarcity of data, especially in the near-threshold region, makes 
this study more difficult. In fact there is only one total cross section and one angular 
distribution available near-threshold at an excess energy of Q = 83 MeV measured by the 
DISTO Collaboration [9]. A theoretical study of the pp — > ppcf) reaction in Ref. [31] was 
made within a similar approach to that of the present study, in the sense that it used a 
relativistic meson-exchange model, considering contributions from the nucleonic and 0p7r 
exchange currents as the dominant contributions. The present study differs from that of 
Ref [31] in that: (1) the angular distribution data from DISTO [9] used in this study 
were reanalyzed [9] and absolute normalization of the corresponding total cross section was 
established, and (2) the pp ISI is included explicitly. 

The relative importance among the possible meson exchange current contributions was 
estimated based on an SU(3) effective Lagrangian, together with various effects described 
in Ref. [31]. Furthermore, the test calculations performed in Ref. [31] showed that the 4>fm- 
exchange current was by far the dominant mesonic current. The com6mec? contribution of all 
other meson-exchange currents to the total cross section is about two orders of magnitude 
smaller than this. Moreover, possible contributions from meson-exchange currents involving 
heavy mesons, in particular, the (f)(f)fi- and 0cj/i- exchange currents were also examined using 
the larger values of the coupling constants calculated from the observed decay of /i — > 0-1-7. 
However, this contribution, as well as the (jxpa- and 0a;(7-exchange currents, also turned out 
to be negligible [31]. Finally, as in the case of cu production, there are neither experimental 
indications of any of the known isospin-1/2 A^* resonances decaying into the N(f) channel, 
nor do there exist enough data for the near-threshold pp ppcj) reaction to fix the relevant 
parameters and judge their validities. Thus, we study the pp — > pp0 reaction considering 
the contributions only from the nucleonic and 0p7r mesonic currents. 

The coupling constant, gf^^Tr, associated with the mesonic current J^^ci ^e extracted 
from the measured branching ratio. Specifically, the coupling constant g^p.,^ = —1.64 is 
determined directly from the measured decay width of ^ p -|- tt [58], where the sign is 



15 



inferred from SU(3) symmetry [31]. We note that the coupling constant g^f^pr^ = —1.64 is 
extracted at different kinematics compared to that of the 00: g^p.^ is determined at = 
and ql = m^, whereas g^pT^ is extracted at = and q^ = m^. Then, the corresponding 
form factor (cf. Eq. (13)) is defined by 

F,p^{ql ql) ^ Fp{ql)FMl) = {^Jf^) (f^) " ^^4) 

In Eq. (24) we again assume the cut-off parameter = Ap — corresponds to those in 
Eq. (13). 

Ahhough we do not have to use the same parameters for the production as those used 
in the pp — > ppur reaction, we use the same value An = 1190 MeV for the cut-off parameter 
in the ^A'^A'^ meson production vertex, because the (pfm exchange current gives the dominant 
contribution to the total cross section, and reproducing the absolute normalization of the 
total cross section is relatively insensitive to the cut-off parameter A^r in the nucleonic 
current compared to the cut-off parameters in the (ppn vertex form factors. Therefore, we 
have three free parameters to be adjusted to reproduce the meson production total cross 
section and angular distribution data [9], namely, gif^NN and in the nucleonic current, 
after replacing a; — > in Eq. (10), and A^ = Ap — A^^ in the form factor of Eq. (24) in the 
(ppn mesonic current, after replacing a; ^ in Eqs. (12) and (14). 

In Figs. 10 and 11 we show the dependence of the calculated angular distributions, 
for g(j)NN = —0.4 and g,j,NN = —1.6, respectively. Results in Fig. 10 imply that as long 
as the value of gt^NN is small, the angular distribution data can be reproduced well within 
experimental error bars, irrespective of the values of up to ~ —4.0. On the other 
hand, results in Fig. 11 show that the larger value, g(j)NN — —1.6 makes the shape of the 
calculated angular distribution sensitive to k^. One can notice that for a certain value of 
K^, the shape of the calculated angular distribution changes from convex to concave. After 
some test calculations, we find the optimum value for this transition is roughly k.^ ~ —2.0. 
Thus, around ~ —2.0 we can expect that there are a large number of possibilities for 
the values oi g^j^^ and A^ which can reproduce the experimentally observed fiat angular 
distribution [9]. 

Next, we fix the value = —2.0 (and Ajv = 1190 MeV), and study the g^NN dependence 
of the angular distribution. Some of the calculated results are shown in Fig. 12. Note that, 
the top- right panel in Fig. 12 has a contribution solely from the mesonic current {g^NN — 0), 
which, neglecting the (f) — uj mixing and the OZl-allowed two step processes [31], may be 
regarded as the limiting case of no ss component in the nucleon wave function, if the value 
of g4>NN is considered as a measure for the ss component. Since we have not included quark 
degrees of freedom explicitly, it is difficult to draw a definite conclusion on the ss component 
in the nucleon wave function. We summarize in Table III the four possible parameter sets 
for g(i,NN,i^<t> and A^ fixed by fitting the DISTO angular distribution data [9]. They all 
reproduce the experimental data [9] reasonably well. This suggests that one needs to study 
additional observables, e.g., the energy dependence of the pp — >■ pp0 total cross section, in 
order to constrain better the parameters of the model. 

Here, it may be interesting to compare the values of g^pNN and K^p obtained in the present 
work with those extracted in Ref. [65] by studying the off-shell time-like nucleon form factors 
using the p(7, e+e")^ reaction. They obtained {g,f,NN-i>^^) — (1.3,7.2). (Note that their 
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definition of re^ is different from that of the present study by a factor ArriN/m^ {rnN,(j,'- 
masses of the nucleon and meson); for comparison, this factor is included in the value of 
here.) 

Next, using the four parameter sets given in Table III, we study the energy depen- 
dence of the pp — > ppcf) total cross section. We show the calculated results in Fig. 13. 

The results exhibit very similar energy dependences for the parameter sets {g<pNN,i^<p) = 
(0.0, 0.0), (-0.4, -0.5) and (-0.4,-4.0), while that for (-2.0,-2.0) shows a different de- 
pendence especially at excess energies in the region Q < 50 MeV. Thus, measuring the 
energy dependence of the total cross section for Q < 50 MeV will help constrain better the 
model parameters, in particular, the magnitude of the coupling constant q^nn- 

VI. SUMMARY AND DISCUSSION 

We have studied the pp ppu / (f) reactions using a relativistic effective Lagrangian at the 
hadronic level, including both the initial and final state pp interactions. For both reactions 
we have made use of the recently measured uj and (p angular distributions in addition to the 
total cross section data to fix the model parameters. 

We have studied the pp — > ppuj reaction considering two possibilities, i.e., the ui meson 
is produced by: (1) the nucleonic and mesonic current contributions, and (2) the nucleonic, 
mesonic and nucleon resonance current contributions. The results show that the energy 
dependence of the total cross section in the range of excess energies Q < 173 MeV, is 
apparently described better by the inclusion of nucleon resonances, which is implemented 
in a way to be consistent with the 7r~p — > un reaction. However, the calculation still 
underestimates the SATURNE data by about a factor of two, where the data points are in 
the range of excess energies Q < 31 MeV. This remains still a problem in understanding the 
reaction mechanism. In this connection, we need more data for exclusive observables in the 
energy region above, but close to Q = 30 MeV because there is no estabhshed method for 
removing the multi-pion background associated with the a;-meson width from the raw data. 
This removal is necessary for extracting the cross sections in the energy region, Q < 30 
MeV, where the effect of the width is very important, and the extraction can be highly 
model dependent. 

In connection with studying resonance contributions to the pp ppuj reaction, we plan 
to investigate the puj invariant mass distributions for this reaction [63]. A measurement 
of the invariant mass distributions for this reaction, should give significant information as 
to whether contributions from resonances arc appreciable or not. Such theoretical studies 
have been made for the pp — >• ppr] [60], and pp — > pAK~^ [66] reactions. Thus, the study 
of the invariant mass distributions may be an alternative method for studying the possible 
a;-meson (and 0-meson) resonance couplings both theoretically and experimentally. 

In addition to the pp ppu reaction, we have studied the pp ppcf) reaction considering 
the contributions solely from the nucleonic and mesonic current contributions. Because of the 
scarcity of data for this reaction in the near threshold energy region, we have obtained four 
parameter sets which can reproduce the (p angular distribution data from DISTO [9] equally 
well. Predictions for the energy dependence of the pp — > ppcf) total cross section indicate 
that a measurement of the cross section close to threshold should be able to constrain better 
the coupling constant g,j,NN (— or ~ —2). 
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Finally, although there exists an enormous interest in vector meson properties in highly 
compUcated many-nucleon environments, e.g., dilcpton production in heavy ion colhsions 
and meson (a;) nuclear bound states, the data for the NN —>■ NNV reaction {V: vector me- 
son) are currently inadequate for understanding the production mechanism of these mesons 
in free space. Thus, more measurements of vector meson production in free space, and espe- 
cially in NN collisions, may be a first step towards understanding the properties of vector 
mesons in such complicated nuclear environments. 
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TABLES 



TABLE I. Model parameters fixed for the pp ppu reaction without the inclusion of nucleon 
resonances. Below, "Bonn" indicates that the same value in the Bonn NN potential B (Table 
A.l) [41] is used. 



Vertex 


Coupling constant 


Cut-off (MeV) 




Nucleonic current: 


[fuiNN = K^guiNN] 






LoNN [lo production] 


guiNN = 9.0 


Aiv = 1190 [See 


Eq. (11).] 




[«a; = -2.0] 






MNN[M = TT, rj, p, Lo, a, ao{= S)] 


Bonn 


Bonn 




Mesonic current: 








Lopn [to production] 


gu,fm = 10.0 




= 1000 [See Eq. (13).] 


pNN 


Bonn 


Bonn 




ttNN [pv-coupling] 


Bonn 


1300 
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TABLE II. Model parameters fixed for the pp ppto reaction with the inclusion of nucleon 
resonances. Below, "Bonn" indicates that the same value in the Bonn NN potential B (Table 
A.l) [41] is used. 



Vcl LtJX 


v^oupiiiig constant 




Put r>fF I'A/ToV'l 

v>UTj-on \^ivie V j 


Nucleonic current: 


\ f n 1 






UJly ly ICC JJiUU-UCUHJllJ 


rt UTAT- O f| 

yujNN — y-u 
— _o 51 






MNN\M = TT 71 n uj CT nn(= 






Bonn 


Mesonic current: 








u/piT [io production] 


inn 

9ujpTr — 10.0 




Ap = ooU, = i4oU [bee iiiq. (^z 


^ AT AT 

pl\ iv 


Bonn 




Bonn 


ttNN [pv-coupling] 


Bonn 




yuu 


Spin 1/2 resonance current: 


qmnnjQmnn* 








mM*+mN \9p,ojNN*, 


l^p,u)gp,ojNN* 


) [See Eqs. (18d) and (18e).] 


S'ii(1535),r = 150 MeV 








MNN[M = TT,ri,p,uj] 


Bonn, but g^iNN = 


9.0 


Bonn, but A.^nn = 900 


[ttNN (pv-coupling)] 








7rArS'n(1535) 


1.25 




900 


r?iVS'ii(1535) 


2.02 




Bonn 


pN Su{l535) 


(0.0, -4.50) [fm] 




Bonn 


ooNSu{1535) 


(-1.04, 3.82) [fm] 




Bonn 


Pii(1710),r = 100 MeV 








]\/f AT AT\ A/f — /T TT n <^ 
IvI iV iV [jkZ — C7, 71, Ty, yy, u/J 


oonn, DUT/ yijjNN — 


Q n 


jDonn, uuTi ^i-T^NN — yuu 


[ttJVTV (pv-coupling)] 










-4.30 




Bonn 


Tr/VPii flTlO) 


1.20 




900 


r7iVPii fl710) 


4.43 




Bonn 


pNPa^ fl710) 


(0.0, 6.70) [fm] 




Bonn 


a;iVPii flTlO) 


fO -1 19) ffml 




Bonn 


bpm 3/2 resonance current: 


(1) 

gMNN,gMNN* 








r (1) / (2) 


= -2.1] 


[See Eqs. (22c) and (22d).] 


1)13(1700), r = 100 MeV 








MiViV[M = 7r,p,uj] 


Bonn, but g^iNN = 


9.0 


Bonn, but A^^nn = 900 


[kNN (pv-coupling)] 








7rAr£>i3(1700) 


0.44 




900 


nNDToi^ 700") 






Bonn 


ojNDi:i{nm) 


3.02 




Bonn 


Pi3(1720),r = 150 MeV 








MNN[M = t:,p,uj\ 


Bonn, but g^jNN = 


9.0 


Bonn, but A^jvAf = 900 


[ttNN (pv-coupling)] 








7rArPi3(1720) 


0.17 




900 


pA^Pi3(1720) 


-3.73 




Bonn 


a;A^Pi3(1720) 


3.94 




Bonn 
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TABLE III. Model parameters for the pp pp(j) reaction, for four possible sets by the </> angular 
distribution, denoted by (a), (b), (c) and (d). Below, "Bonn" indicates that the same value in the 
Bonn NN potential B (Table A.l) [41] is used. 



Vertex 


Coupling constant 




Cut-off (MeV) 


Nucleonic current: 


[frj>NN = fi(j)9<j)NN] 






(f)NN [(f) production] 


(a) g4,NN = 0.0 [k^ = 0] 




Ajv = 1190 [See Eq. (11). 




(b) gci>NN = -0.4 [k^ = 


-0.5] 


Ajv = 1190 




(c) g4,NN = -0.4 [k^ = ■ 


-4.0] 


Ajv = 1190 




(d) g4,NN = -2.0 [kci> = 


-2.0] 


Ajv = 1190 


MNN[M = -K,7],p,LV,a, ao(= 6)] 


Bonn 




Bonn 


Mesonic current: 






A^ = Ap = A^ 


(a) 0p7r [(p production] 


gcppTT = -1-64 




A^ = 1930 [See Eq. (24).] 


(b) 4>pTT [(j) production] 


g4>fm = -1-64 




A,/, = 2100 


(c) (ppiT [(/) production] 


g<f>pn = -1.64 




A^ = 1915 


(d) 4>pn [(f) production] 


g(i>pn = -1.64 




A,^ = 2200 


pNN 


Bonn 




Bonn 


ttNN [pv-coupling] 


Bonn 




1300 
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FIGURES 





FIG. 1. Decomposition of the reaction amplitude for the NN — ^ NNV {V = u, 0) reaction in 
the present work. Tmn denotes the meson-nucleon T-matrix. ISI and FSI stand for the initial and 
final state NN interactions, respectively. 



25 




FIG. 2. Vector meson {u or (f>) production currents, J'^, included in the present study: a) 

nucleonic and resonance currents, b) mesonic current. V = uj oi (j) and M = TT,r}, p,uj,a,ao{= S). 
In the intermediate states of diagram a) , negative-energy propagations for nucleon and resonances 
are included. 
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-1 -0.5 0.5 -1 -0.5 0.5 1 

cos(e) cos(0) 

FIG. 3. dependence of the uj angular distribution at excess energy Q = 173 MeV, without 
the inclusion of nucleon resonances. The (dashed, dot-dashed, solid) lines show the (mesonic, 
nucleonic, total) contributions, respectively. The dots denote data from COSY-TOF [38]. The 
cut-off parameter Ajv in the form factor at the uj production vertex uNN is fitted to the total 
cross section of 30.8/u6 at excess energy Q = 173 MeV [38] for each panel. 
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-1 -0.5 0.5 -1 -0.5 0.5 



COS(0) COS(0) 

FIG. 4. g^jNN dependence of the lo angular distribution at excess energy Q = 173 MeV. The 
left panel is one of the reasonable fits achieved with g^jMN = 9.0, while the right panel is the result 
obtained with the value g^iNN = 17.37, which is approximately the value used in the Bonn NN 
potential model [41]. Both calculations use the value, = —2.0. Also, see the caption of Fig. 3. 
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10^ 
10^ 
10^ 
10° 
10-^ 

!q 10"^ 

D 10' 

10^ 
10° 
10" 

10"" 



I ^ 

FSI+width 










1 1 1 1 MM 


1 1 1 1 Mil 


h- H-4 



Phase space + FSI 
(Bonn) 




- no width 

— width 



(normalized at Q=92 and 173 MeV) 



— I — I I I 1 1 1 1 1 — I — I I I 1 1 1 1 

FSI (no width) 



rrrri r- 




Nucleonic 
Mesonic 
Total 



ISI+FSI+width 




Hibou et al. 



10^ 



lO'' 



10 10 10-^ 10 10 10 10 
Q [MeV] Q [MeV] 

FIG. 5. Energy dependence of the total cross section for the pp — > ppco reaction without the 
inclusion of resonances. "ISI", "FSI" and "width" stand for the pp initial state interaction, pp 
final state interaction, and effects of the u width, respectively. The result with "ISI+FSI+width" 
(the bottom-right panel) should be compared with the data, where those panels without any of 
the legends, "ISI" , "FSI" and "width" , imply that the corresponding effect is switched off from the 
full calculation (the bottom-right panel). Data are from SATURNE [37] (dots), COSY-TOF [38] 
(diamonds), Ref. [36] (circles) and DISTO [9] (a filled square), respectively. "Hibou et al." (the 
bottom-right panel) stands for the result used in the analysis in Ref. [37] . In the bottom-left panel 
indicated by "Phase space + FSI", the calculated energy dependences are normalized to the total 
cross section data from COSY-TOF [38] at cither Q = 92 or 173 MeV, where "(e* • J) = constant" 
is used in Eq. (1), and thus, effects of the FSI generated by the Bonn NN potential model [41] 
enter. 
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10 



71 + p ^ CO + n 



Nucleonic 

Mesonic 

Nuc+Mec 

— N*jj sum 

N*j3 sum 

Resonance sum 

Total 



A^=1100 MeV, A =850 MeV, n =2 





1.70 




1.75 



1.80 



1.95 



2.00 



2.05 



1.85 1.90 
W (GeV) 

FIG. 6. Energy dependence of the total cross section for the tt" + p ^ un reaction obtained 
with the preferred model parameter set. Data are from Ref. [53]. Note that at an excess energy of 
Q = 173 MeV in the pp ppu reaction, the maximum center-of-mass energy W for this reaction 
reaches c± 1.9 GeV. 
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cos(e) 
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fA AJ=(0.85,1.45) GeV 



• COSY-TOF 

Nucleonic 

Mesonic 

Resonances 
Total 




-1.0 



-0.5 



0.5 



1.0 



173 



0.0 

cos(e) 

FIG. 7. Calculated uj angular distribution for the pp ppiv reaction at excess energy Q 
MeV. The cut-off parameter Ajv in the loNN form factor, Fj\[{p'^) = A^/[A^ + (p^ — m^)] 
[Eq. (11)], is adjusted to the Tr^p ujn reaction Ajv = 1100 MeV (Fig. 6), and the cou- 
pling constant gaNPn is fitted to reproduce the total cross section of 30.8/x6 at Q = 173 MeV. 
The cut-off parameters and Hp indicated in each panel enter at the ujpiT vertex form factor: 
F.p.{ql,q'p) ^ Fp{ql)xF^{ql) = [A2/(A2 - q^^ >c[{Kl - ml)/{Kl q^)] [Eq. (23)], withn, = 2, 
Ap = 850 MeV and Att = 1450 MeV, respectively. Also, see the caption of Fig. 3. 
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N"^ sum 

Total 



I I I I I I I I L 
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Q [MeV] 



10 
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FIG. 8. Calculated energy dependence of the total cross section for the pp ppu reaction 
with the inclusion of nucleon resonances, Sn (1535), Pn (1710), Di3 (1700) and Pi3(1720). Also, see 
the caption of Fig. 5. 
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FIG. 9. Decomposition of resonance contributions for the energy dependence of the — > ppuj 
total cross section. Also, see the caption of Fig. 5. 
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K^=0.0 

A^=2150 MeV 
g^,,=-0.4,AK=1190 MeV 




K^=-4.0 
A.=1915 MeV 




K^=-0.5 
A.=2100 MeV 



-1 



-0.5 





cos(e) 



0.5 



-1 



-0.5 





cos(e) 



0.5 



FIG. 10. dependence of the (j) angular distribution at Q = 83 MeV. The (dashed, dot-dashed, 
sohd) hues show the (mesonic, nucleonic, total) contributions, respectively. The <^NN coupling 
constant and the cut off parameter A^r associated with the (pNN meson production vertex are 
fixed at g^NN = —0.4 and Ajv = 1190 MeV, and the cut-off parameter = Ap = A^^ in the (pp'K 
meson production vertex is fitted to reproduce the total cross section of 190n6 at Q = 83 MeV. 
The dots are the data from DISTO [9] (a). 
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FIG. 11. Same as Fig. 10, but q^nn = —1.6. 
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FIG. 12. g(f)NN dependence of the (p angular distribution with the fixed value n^p = —2.0, and 
An = 1190 MeV. Also, see the caption of Fig. 10. 
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FIG. 13. Energy dependence of the total cross section for the pp ppcf) reaction, calculated 
using the four parameter sets, which can reproduce the (j) angular distribution data [9]. A systematic 
error in the data point from DISTO [9] is not included. 
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